Conditional generation of arbitrary single-mode quantum states of light 

by repeated photon subtractions 
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We propose a scheme for the conditional generation of arbitrary finite superpositions of (squeezed) 
Fock states in a single mode of a traveling optical field. The suggested setup requires only a source 
of squeezed states, beam splitters, strong coherent beams, and photodetectors with single-photon 
sensitivity. The method does not require photodetectors with a high efficiency nor with a single- 
photon resolution. 
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I. INTRODUCTION 

During recent years it has been widely recognized that 
non-classical states of light represent a valuable resource 
for numerous applications ranging from ultra-high pre- 
cision measurements pl 0, to quantum lithography 
0, Ü| and quantum information processing |(| . It is of- 
ten desirable to generate nonclassical states in traveling 
optical modes, as opposed to the cavity QED experi- 
ments where the generated state is confined in a cav- 
ity and can be probed only indirectly. Many ingenious 
schemes have been proposed and experimentally dcmon- 
strated to generate the single-photon states Q, 13 and 
various multiphoton entangled states such as the GHZ 
states 0, eluster states |Ï3, and the so-called NOON 

states (nimEiQm. 

Considerable attention has been also devoted to the 
preparation of arbitrary single-mode states 0, 0, 0, 
l·LS| and, in particular, the Schròdinger cat-like superpo- 
sitions of coherent states (2(3, [^J which can represent 
a valuable resource for quantum information processing 
[22L l23| . Experimental generation of arbitrary superposi- 
tions of vacuum and single-photon states has been accom- 
plished using a parametric down-conversion with input 
signal mode prepared in a coherent state |24|, employ- 
ing the quantum scissors scheme pM |26| or conditioning 
on homodyne measurements on o ne part of a non-local 
single photon in two spatial modes [27|. It is, however, 
very difficult to extend these experiments to superpo- 
sitions involving also two or more photons. The known 
schemes for conditional generation of arbitrary superposi- 
tions of Fock states require single-photon sources and/or 
highly efficient detectors with single photon resolution 
which represents a formidable experimental challenge. 

In this paper, we propose a novel state preparation 
scheme, which does not require single-photon sources and 
can operate with high fidelity even with low-efficiency 
detectors that only distinguish the presence or absence 
of photons. Our scheme is inspired by the proposal of 
Dakna et al. |16| who showed that an arbitrary single- 
mode state can be engineered starting from vacuum by 
applying a sequence of displacements and single-photon 
additions. Our crucial observation is that if the initial 



state is a squeezed vacuum, then the single-photon ad- 
dition can be replaced with single-photon subtraction 
[28L I29I ]. which is much more practicable. Indeed, a 
single-photon subtraction can be achieved by diverting 
a tiny faction of the beam with a beam splitter towards 
a photodetector, so that a click means that a photon has 
been subtracted from the beam (this process becomes 
exact for a transmittance tending to one). In fact, the 
single-photon subtraction from a squeezed vacuum has 
already been experimentally demonstrated |30| . which 
provides a strong evidence for the practical feasibility 
of our scheme. We note that the photon subtraction is 
an extremely useful tool that allows to generate states 
suitable for the tests of Bell inequality violation with 
balanced homodyning |3lL |32j , can be used to improve 
the performance of dense coding [33j, and forms a cru- 
cial element of the entanglement distillation schemes for 
continuous variables |34j . 

The present paper is organized as follows. In Sec- 
tion II, we explain the mechanism of state generation on 
the simplest non-trivial example of preparing a squeezed 
superposition of vacuum and single-photon states. Our 
setup then consists of two displacements and one condi- 
tional photon subtraction. We present the details of the 
calculation of the Wigner function of the generated state 
for a realistic setup involving imperfect photon subtrac- 
tion (obtained with imperfect detectors and beam split- 
ters with a non-unity transmittance). In order to eval- 
uate the performance of the scheme, we investigate the 
achieved fidelity and the preparation probability for var- 
ious target states. In Section III, we extend the scheme 
to the generation of an arbitrary single-mode state and 
show how to calculate the displacements that need to 
be applied during the state preparation. As an illustra- 
tion, we consider the generation of several states which 
are squeezed superpositions of vacuum, single-photon, 
and two-photon Fock states. Finally, the conclusions are 
drawn in Section IV. 

II. GENERATION OF A SQUEEZED 
SUPERPOSITION OF |0) AND |1) 

In this section, we introduce our setup for the gener- 
ation of an arbitrary squeezed superposition of vacuum 



2 



and single-photon state, which consists of two displace- 
ments with a photon subtraction in between, as schemat- 
ically sketched in Fig. 1. This setup represents a bàsic 
building block of our universal scheme: as shown in Sec- 
tion III, any squeezed superposition of the first N + 1 
Fock states can be generated from a single-mode squeezed 
vacuum by a displacement followed by a sequence of N 
photon subtractions and displacements. 



A. Pure-state description 

We first provide a simplified pure-state description of 
the setup, assuming perfect detectors with single-photon 
resolution, which will give us an insight into the mecha- 
nism of state generation. We will show that conditionally 
on observing a click of the photodetector PD, the setup 
produces a squeezed superposition of vacuum and single- 
photon states, 



PD 



/target — S(a)[co|0)+ci|l>], 



(1) 



where S(s) — exp[s(a t2 — a 2 )/2] denotes the squeezing 
operator, and a (cr) is the annihilation (creation) oper- 
ator. 

Our state engineering procedure starts with a single- 
mode squeezed vacuum state, which is generated in an 
optical parametric amplifier (OPA), 



S(Sin)|0> 



f 



Vcosh( Sin ) ^ 2«n! 

where s; n denotes the initial squeezing constant. The 
single-mode squeezed vacuum passes through three 
highly transmitting beam splitters, which realize a se- 
quence of a displacement followed by a single-photon sub- 
traction and another displacement. The state is displaced 
by combining it on a highly unbalanced beam splittcr 
BS_d with transmittance To > 99% with a strong coher- 
ent state \a/ru), where rn_ = VÏÏd and Ro = 1 — Td 
is the reflectance of BSd [35j. In the limit To — > 1, the 
output beam is displaced by the amount a. This method 
has been used, e.g., in the continuous- variable quantum 
teleportation experiments j3(|. For the sake of simplic- 
ity, we shall assume that To = 1 and the displacement 
operation is exact. The conditional single-photon sub- 
traction requires a highly unbalanced beam splitter BS 
with transmittance T, followed by a photodetector PD 
placed on the auxiliary output port. A successful photon 
subtraction is heralded by a click of the detector. In the 
limit T — ► 1, the most probable event leading to a click 
of the detector is that exactly a single photon has been 
refiected from the beam splitter. The probability of rc- 
moving two or more photons is smaller by a factor oïl—T 
and becomes totally negligible in the limit T — > 1. The 
conditional single-photon subtraction can be described 
by the non-unitary operator 
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FIG. 1: Proposed experimental setup for generating \tp) = 
S'(s)[co|0) + ci 1 1)] . An optical parametric amplifier (OPA) 
generates a single-mode squeezed vacuum state which then 
propagates through three highly unbalanced beam splitters 
BSd, BS, and BSd, which realize a sequence of two displace- 
ments interspersed with one conditional photon subtraction. 
Successful state preparation is heralded by a click of the pho- 
todetector PD. 



where n — a^a is the photon- number operator, while t = 
yT and r = \/l — T denote the amplitude transmittance 
and reflectance of BS, respectively. 

The input-output transformation corresponding to the 
sequence of operations in Fig. I reads 



/out — D{a 2 )XD(ax)S(s in )\0), 



(4) 



where D(a) = exp(aa^ — a*a) is the displacement oper- 
ator. We will show later on how the displacements a.\ 
and «2 depend on the target state 0J, as well as how the 
initial squeezing Si n depends on the target squeezing s 
for a given transmittance T < 1. But, to make it simple, 
let us assume first that T = 1 and ai = —a 2 = a. Then, 
using D(ay aD(a) = a + a, the conditionally generated 
state can be written as 



)out = (a + a)S(s in )\0}. 



(5) 



Taking into account that a and a' transform under the 
squeezing operation according to 



(s)aS(s) — a cosh(s) + sinh(s) , 
S , ' i '(s)a^S , (s) = a) cosh(s) + asinh(s), 



(6) 



we obtain 



|-0)out = S , (si n )[acosh(s in ) + a f sinh(s in ) + a]|0) 

= S(s in )(smh( Sin )|i) + a|0)). (7) 

We can see that by setting a = (cq/ci) sinh(si„), we ob- 
tain the target state JTJ. This simple analysis illustrates 
the principle of operation of the scheme shown in Fig. 
1. However, the limit T = l is unphysical, because the 
probability of successful state generation vanishes when 
T — > 1. Let us now take into account T < 1. 

In order to simplify the expression we first rewrite 
all displacement operators in a normally-ordered form, 



D(a) 



*l72 f 



a , and we obtain 



X = t n ra, 



(3) 



)out oc e Q2a V a ^í™ae Qia V a '" 



S(Sin)|0>. (8) 
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Next, we propagate the operator t n to the right by using 
the relations 



t n e a * a = e a * a/t t n , t n e aa = e taa t n . (9) 



After these algebraic manipulat ions we obtain 

|^>out oc e Q2at ae tQlQt e- [ ^ +Q i /tla í™S*(s in )|0). (10) 
Note that we have also moved to the right the operator 



e "í" and used the fact that e 



The combined action of the operators t n S(si n ) on vac- 
uum produces a single-mode squeezed vacuum statc 
just as without applying í n , only with a lower squeezing 
constant s satisfying tanh(s) = í 2 tanh(s in ). Thus, we 
can write 



t n S(s in )\0) ocS(a)|0>. 



(11) 



Finally, we move the operator e Q2<1 to the right, using 
the formula e" 2a+ a = (a — a 2 ) e Q2at , which results in 



|^) ou tcx (a-a 2 )e Sa 'e^ a S(s)\0), 



(12) 



where S = «2 + toi and 7 = 0:2 + cti/í. With the help of 
Eq. JïïJ), we can write, 

S(s)\0) OC 5( s ) e [-5cosh( S )- 7 * S inh( S )]a+| ^ ^3) 



„<5a 1 ' 



-7 a 



which is a state with a generally non-zero coherent dis- 
placement. This displacement vanishes provided that ot\ 
and a 2 satisfy 

(cü2 + toi) cosh(s) = («2 + &i/t) sinh(s). (14) 

If the condition l|14(l holds, then the output reads 



> out oc (a- a 2 ) S(s)\0) 

oc S , (s)[a' 1 ' sinh(s) — a 2 ] 



(15) 



where we have moved the squeezing operator to the left 
using Eq. JïïJ). The desired squeezed superposition of the 
first two Fock states can then be obtained by choosing 



«2 



ço 

Cl 



sinh(s) 



(16) 



[tanh 2 (s) - t 2 ]a 2 + (t 2 - 1) tanh(s)c4 

«i — t - q , (17) 

í 4 -tanh 2 (s) V ; 

where the displacement a\ is determined from the condi- 
tion (|14[) . Note that we may assume that the coefficient 
ci of the Fock state |1) is non-zero; otherwise, no photon 
subtraction is needed to generate the target state. 

It should be emphasized that the generated state with 
our scheme is a squeezed superposition of Fock states. 
In principle, this squeezing could be removed by sending 
the conditionally generated state (|15f) through a squeezer 
S{—s). However, in many cases, the squeezing may not 
be an obstacle or may even represent an advantage, such 
as in the generation of Schròdinger cat states | a) — \ — a) 
which can be for small \a\ very well approximated by a 
squeezed single-photon state S(s) |1) [2lll29| . 



B. Realistic model 

We shall now present a more realistic description of the 
proposed scheme taking into account that the photode- 
tectors exhibit only single-photon sensitivity, but cannot 
resolve the number of photons in the mode, and have 
a detection efficiency 77 < 1. Such detectors have two 
outeomes, either a click or a no-click. We model this 
detector as a sequence of a beam splitter with transmit- 
tance 77 followed by an idealized detector which performs 
a projection onto the vacuum and the rest of the Hilbcrt 
space, n = |0)(0| (no click), lli = 1 - |0)(0| (a click). 

Similarly as in Ref . |3l| , it is convenient to work in the 
phase-space representation and consider the transforma- 
tion of Wigner functions. The setup in Fig. 1 involves 
two modes, the principal mode A and an auxiliary mode 
B. The Gaussian Wigner function of the initial state of 
mode A after the first displacement reads 



W G (r A ;T A ,d A ) 



v / detTI _ 



(r A -d A ) r A (r A -d A ) 



(18) 



where r A = (x A ,p A ) T is the vector of quadratures of 
mode A and d A = z\ = -\/2(5toi, Sai) T is the dis- 
placement. The matrix T A is the inverse of the covari- 
ance matrix 7^4. Initially, the mode A is in a squeezed 
vacuum state and the covariance matrix is diagonal, 
7.4 = diag(e 2s '=,e- 2s '=). 

In a second step, the modes A and B are mixed on 
an unbalanced beam splitter BS and then mode B subse- 
quently passes through a (virtual) beam splitter of trans- 
mittance r\ which models the imperfect detection with ef- 
ficiency 77. This transformation is a Gaussian completely 
positive (CP) map M, and the resulting state of modes A 
and B is still a Gaussian state with the Wigner function 



WabÍtab) = Vdet , rAfl e-( f ^-43) r r, B ( w -^} | 

(19) 

where r A B — (x a ,p a ,xb,Pb) T ■ The vector of the first 
moments d A s = (d A ,dB) T and the covariance matrix 
7ab = can be expressed in terms of the initial pa- 
rameters of mode A before the mixing on BS (i.e., z\ and 
~y A ) as follows: 



d A i 



d A 

<ÍB 



1AB = S{ 1A ®I B )S T + G, 



(20) 



where S — S v SbSi S n = I A © ^/t]Ib and G = A 
(1 — tÍ)Ib model the inefficient photodetector, and Sbs is 
a symplectic matrix which describes the coupling of the 
modes A and B on an unbalanced beam splitter BS, 



Sbs 



t r 

t r 

-r t 

-r t 



(21) 
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FIG. 2: (a) Fidelity between the generated state and the 
target state and (b) probability of successful generation as 
a function of the squeezing Si n for the four target states 12911 
(sòlid line), lEÏÏJl (dashed line), (ISH (dot-dashed line) and lO 
(dotted line), with T = 0.95 and r) = 0.25. 



After the photon subtraction, the density matrix pA,out 
of mode A conditioned on a click of the photodetector PD 
measuring the auxiliary mode B becomes 



(22) 



where Trs denotes a partial trace over mode B, and pab 
is the two-mode density matrix of the Gaussian state 
characterized by the Wigner function (|19fl . Then, af- 
ter the second displacement of z 2 = V2(3ia 2 ,^sa 2 ) T , the 
Wigner function of mode A can be written as a linear 
combination of two Gaussian functions i|18|) , namely 

W(r) P = CiW G (r; Ti,d x ) + C 2 W G (r; T 2 ,d 2 ), (23) 

where P is the probability of successful generation of 
the target state. The expression (|2l5|) can be derived by 
rewriting Eq. (1221 in the Wigner representation. One 
uses the fact that the Wigner function of the POVM el- 
ement Hi,B is a diffcrence of two Gaussian functions, 



FIG. 3: (a) Fidelity between the generated state and the 
target state and (b) probability of successful generation as 
a function of T for the four target states 129I I- 132II . The 
curves are plotted considering the optimal squeezing Si n for 
each state, namely 0.50 dB for state l|29|l (sòlid line), 1.66 dB 
for state 13011 (dashed line), 0.85 dB for state 1311 1 (dot-dashed 
line), and 0.36 dB for state ^ (dotted line). The curves are 
plotted for rj = 0.25. 



Similarly, the formulas for the parameters of the second 
tcrm read 



T, = T A -crT 



3 O , 



d, 
C-i 



cIa + T 2 1 aT B 1 dB + z 2 

<1<''ÍI ,w; : . ,, 



V det(r 2 )det(T B ) 
where Tb = ^b + I and 



cxp [-d T B Md B ] , (27) 



M = ÏRÏ 



B± B 



(28) 



Since all the Wigner functions in Eq. Il2.''>ll are normalized, 
the probability of a successful state generation can be 
calculated simply as the sum P = C\ + C 2 . 



1 

2^ 



1 



-X —p 



and that the trace of the product of two operators can 
be evaluated by integrating the product of their Wigner 
representations over the phase space. 

To dcfinc the matrices and vectors appearing in 
Eq. (|23|l . we first divide the matrix T ab — Iab 
four sub-matrices with respect to the A\B splitting, 



^24) C. Examples 

In order to illustrate our method, let us consider the 
preparation of the following four squeezed superpositions 
of |0) and |1) states, 



^ab = 



Ta cr 
t y 



(j 



B 



(25) 



The correlation matrix Ti and the displacement d\ ap- 
pearing in the first term on the right-hand side of Eq. Ij23(l 
are given by 



di 



d A + z 2 , 
1. 





= S(s)\l), 






= S(s)^=(\0) + 






= S(s) ï_(3|0) 

Vïo 






= S(s) ÍL(|0)-f 

Vïo 


-3|1)). 



(29) 
(30) 

(31) 
(32) 



The fidelity of the generated state for the target states 
(|2~9Tl ll32l is plotted in Fig. 2(a) as a function of squeez- 
ing. We can see that the conditionally prepared states are 
close to the dcsircd states and their optimum fidelities are 
(26) reached for a low initial squeezing (below 2 dB), which is 
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FIG. 4: (a) Fidelity between the generated state and the tar- 
get state and (b) probability of successful generation as a func- 
tion of r/ for the four target states J29I I- J32II . The curves are 
plotted considering the optimal squeezing Si n for each state, 
namely 0.50 dB for state l|29|l (sòlid line), 1.66 dB for state 
l!50t (dashed line), 0.85 dB for state J^D (dot-dashed line), 
and 0.36 dB for state (dotted line). The curves are plot- 
ted for T = 0.95. 

experimentally accessible. Although it is hardly visible 
in Fig 2(a), there is typically a non-zero optimal value 
of the squeezing, giving the highest fidelity. As shown 
in Fig. 2(b) the increase of the squeezing improves the 
probability of successful generation of the target state. 
A comparison of Fig. 2(a) with Fig. 2(b) reveals a clear 
trade-off between the achievable fidelity and the prepa- 
ration probability. 

Figure 3(a) shows the dependence of the fidelity on the 
beam splittcr transmittancc T, considering the optimal 
squeezing for each of the states. (Note that for state (|29ü . 
we could not find numerically the optimum squeezing, so 
we arbitrarily chose s- m — 0.50 dB as an optimal value in 
other to keep a reasonable generation probability.) We 
see that as T approaches unity, the fidelity gets arbitrar- 
ily close to unity, while the probability of successful state 
generation P oc (Í—T) r\ decreases, as shown in Fig. 3(b). 
The value T = 0.95 used in Fig. 2 seems to be a reason- 
able compromise between the success rate (P » 10~ 3 or 
P w 10~ 4 depending on the target state) and the fidelity, 
F > 0.95%. 

We also have studied the dependence of the fidelity 
on the detection efnciency r\. The numerical results are 
shown in Fig. 4(a), where we can see that the scheme is 
very robust in the sense that the fidelity almost does not 
depend on rj. Fidelities above 95% could be reached even 
with r\ of the order of a few percent if T is high cnough. 
This is in agreement with the findings of Ref. [31j . How- 
ever, a low r\ reduces the preparation probability, as 
shown in Fig. 4(b). 

III. ARBITRARY SINGLE-MODE SQUEEZED 
STATE 

In the preceding section, we have demonstrated that 
the combination of two displacements and a photon sub- 
traction allows us to build any squeezed superposition 



of |0) and |1) states. In this section, we shall generalize 
this procedure to any squeezed superposition of the first 
N + 1 Fock states, 

N 

IV'Harget = S(s) }^ C n \n) , (33) 
n=0 

and show that it can be prepared from a squeezed vac- 
uum by applying a sequence of iV+1 displacements inter- 
spersed with N photon subtractions, as shown in Fig. 5. 



A. Pure-state description 

As in the preceding section, we first provide a simpli- 
fied pure-state description of the setup, assuming per- 
fect detectors with single-photon resolution. This will 
allow us to determine the dependence of the coherent dis- 
placements a>j on the target state l|33l) . Gcncralizing the 
procedure presented in the preceding section, the input- 
output transformation corresponding to the sequence of 
operations in Fig. 5 reads 

l^out = D(a N+1 )XD(a N )X . . . D(a 2 )XD(a 1 )S(s ia )\0). 

In order to simplify this expression, we first rewrite all 
displacement operators in a normally-ordered form and 
then move all the operators t n to the right using the rela- 
tions (JÜJ . This results in the substitution ctj — ► otj 
and a* — > a*í J_Ar_1 in the exponents. Next, we prop- 

_fj — N— 1 n ,* n 

agate all the exponential operators e i to the 

right, 

|V>)out oc e aN+ia \é aNa \ . . . ae t " Qiat e- 7 * a í Ar ™5(s in )|0), 

(35) 

where 7 = J2j=i OijP~ N ~ l . The combined action of the 
operators t Nn S(si n ) on the vacuum produces a single- 
mode squeezed vacuum state, t Nn S(s[ n )\0) oc 5(s)|0), 
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FIG. 5: Proposed experimental setup. An optical parametric 
amplifier (OPA) generates a single-mode squeezed vacuum 
state which then propagates through 2iV + 1 highly unbal- 
anced beam splitters BSd and BS, which realize a sequence 
of N + 1 displacements interspersed with iV conditional pho- 
ton subtractions. Successful state preparation is heralded by 
clicks of all photodetectors PD^. 
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where tanh(s) = t 2N tanh(si n ). After some algebraic ma- 
nipulations we get 



JV 



) out oc SO) JJ(acosh(s) + a f sinh(s) - /3j)\0), (36) 



where 



N+l 

& = E ^t N+i - k . 



(37) 



The formula (|3*üjl is vàlid provided that the overall dis- 
placement is zero. This corresponds to the constraint 



N+l 

cosh(s) ajt 

3=1 



JV+1 

x - '-í = sinhO) V a!*»'-*- 1 



3 = 1 



(38) 



which generalizes condition 1)14(1, 

We now prové that an arbitrary superposition of the 



first N + l Fock states ^n=o c "l n ) can be expressed as 



Uf=M-P 3 )\0} = ELo h kA k \0), where A = acosh( S ) + 
sinh(s) and h k are the coefficients of the characteristic 
polynomial whose roots are From the condition 



JV 



JV 



(39) 



k=0 



n=0 



we can immediately determine the coefficients hpf and 
hpf-x, because only the term A N gives rise to a' N and, 
similarly, only the expansion of A 1 ^^ 1 contains a 1 1 . 
We thus get equations 



h N = 



CArsinh N (s) 

7m : 



ÍN-l 



Cjv-isinh 1 N (s) 



(40) 

Once we know and h^-i, we insert them back in 
Eq. JH, and, from Ek=o h kA k \0) = E^o c » " 
{hN-iA 1 * 1 ^ 1 +hjyA N )\0), we determine /ijv-2 and tiN-3- 
By repeating this procedure, we can find all coefficients 
hj. This proves that the condition (|39|l can be always 
met for any nonzero squeezing, hence our method is in- 
deed universal and allows us to generate arbitrary super- 
positions. After finding the hj's, the coefficients (3/s are 
calculated as the roots of the characteristic polynomial 
J2k=o ^ k P k , an( i' finally, the N + l displacements ctj's 
are determined by solving the system of N + 1 linear 
equations (|3T|l and (J3SJ. 



B. Realistic model 



conditioned on a click of photodetector measuring the 
auxiliary mode B is related to pk—i,A as follows, 

p k , A =D k+1 Tr B [M( P k-i,A®\0) B (0\)(lA®'n.i,B)]Dl +1 , 

(41) 

where D k +\ = D(a k +i) is a displacement operator and 
M. denotes the Gaussian CP map (|20|) that describes 
mixing of the modes A and B on BS and accounts for 
imperfect detection. Since each step (|4*T)l gives rise to a 
linear combination of two Gaussian states from a Gaus- 
sian state, the Wigner function of the state pu,A can be 
written as a linear combination of 2 k Gaussian functions, 



2" 

W k {r)P k = J2C j , k W G (r;T jtk ,d j , k ), 

3=1 



(42) 



where P k is the probability of success of the first k pho- 
ton subtractions. The correlation matrices Fj tk and dis- 
placements dj.k after k photon subtractions and k + 1 
displacements can be expressed in terms of Tj, k -i and 

Similarly as in Section IIB, we first define the real 
displacement vector z k = ^/2(%ta kl ( <sa.k) T and the two- 
mode covariance matrix and vector of mean vàlues after 
the action of the CP map Ai , 



Vj,k,A ] = s 
V j,k,B 



2 j,k 

Q 



Jj,k,AB 



>Ib)S t 



G. 



(43) 



We also decompose the inverse matrix Tj t k,AB = 7j k ab 
similarly as in Eq. 125|) . 



r 



j,k,AB 



"fj,k,A a j.k 



T 



(44) 



The J-th term in Eq. I|42|) gives rise to two new terms. 
The (2j — l)-th term is parametrized by 



2j— l,fc — 1 j,k-l,A — a j,k-l 1 



j,k-l,B v j,k-l' 



d<2j-l,k 
C 2 j-l,k 



Vj,k-1,A + z k+l, 

Cj.k-i- 



(45) 



Similarly, the formulas for the 2j-th term read 

Ú2j,k = Vj,k-1,A + ^2j,k a Ó,k-l^J t k-l,B V 3,k-l,B + Zk+1, 
C2j,k = —2Cj.k-l> 



úet(T jt k-i, A B) 



det(r 2 j,fc) det(Tj,fc-i : B) 
x exp [-vT k _ líB Mvj,k-i,B] , 



(46) 



We shall now present a more realistic description of 
the proposed scheme, which takes into account realistic 
photodetectors. After k-th photon subtraction and k + 
1-th displacement, the density matrix pk,A of mode A 



where Tj,fc-i,.B = T jt k-i,B + I and 



M 



j,k-l,B L j,k-1,B 



— 1 T 1 fy\ — 1 

~ 1 j,k-l,B a j,k-l L 2j,k a 3,k-l 1 j,fc-l 5 
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Using these formulas repeatedly starting from the initial 
(k = 0) Gaussian state l|18fl . one obtains after N itera- 
tions the Wigner function of the conditionally generated 
state. The probability of state preparation can be calcu- 

2 N 

lated simply as the sum of Cj n, P = 2J 7=1 Cj n- 



C. Examples 

We shall now consider, as an illustration, the genera- 
tion of squeezed superpositions of |0), |1), and |2). These 
states, namely, 



= S(s) 



1 



y/ï- 



|c | 2 + |ci| 2 



(co|0) + Cl |l) + |2)). (47) 



can be prepared with two photon subtractions. Here, 
we assume that the coeficient c 2 of the Fock state |2) is 
non-zero (we arbitrarily take it equal to one). Otherwise, 
only one (or zero) photon subtraction would be needed to 
generate the target state. In the case of two photon sub- 
tractions interspersed with three displaccments, Eq. Ij36(l 
reduces to 

|V>)out oc 5(s)[(sinh(s)cosh(s)+ / 9 1 /3 2 )|0) (48) 
- (fil +/3 2 )sinh(s)|l) + V2sinh 2 (.s)|2)]. 

This state matches the target state (|47|) if 



01,2 = 



-B ± y/B 2 - 4C 



(49) 



where 



B = v / 2sinh(s)ci, 

C = \/2 sinh 2 (s)c — sinh(s) cosh(s). 

Equations (|37() and l|38|) allow us to calculate the dis- 
placements needed to generate this state. Assuming for 
simplicity that cq, c\ and s are chosen such that fi\ and 
(3 2 are both real, we obtain, 

a 2 = (0i-a 3 )/t, (50) 
tanh(s) (as + a 2 /t) — (a 3 + ta 2 ) 



«i 



t 2 -tanh(s)/í 2 



In order to illustrate our method, let us consider the 
following four squeezed superpositions of the Fock states 



10} , |1>, and |2>: 

m 


= S(s)\2), 


H2», 


(51) 
(52) 




= S(s)-L(\0) - 


H2», 


(53) 




= ^(*)^(|0) - 


Hl> + |2», 


(54) 



1 

0.95 
0.9 

F 

0.85 
0.8 
0.75 



10 



10 



p 

10" 



10" 



2 4 6 

Squeezing(dB) 



10 




2 4 6 
Squeezing(dB) 



FIG. 6: (a) Fidelity between the generated state and the 
target state and (b) probability of successful generation as a 
function of the squeezing Si n for the four target states 151H 
(sòlid line), 1521 (dashed line), 15311 (dot-dashed line), and 
l|54jl (dotted line), with T = 0.95 and 77 = 0.25. 



a) 




0.9 
0.8 
0.7 
0.6 



0.5 u 

0.8 0.85 0.9 j 0.95 1 
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10" 



10" 



10 
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b) 










4 



0.8 0.85 0.9 T 0.95 1 



FIG. 7: (a) Fidelity between the generated state and the 
target state and (b) probability of successful generation as 
a function of T for the four target states 1511 - 1541 . The 
curves are plotted considering the optimal squeezing Si n for 
each state, namely 3.54 dB for state l|51|l (sòlid line), 4.02 dB 
for state 15211 (dashed line), 2.43 dB for state 1531 1 (dot-dashed 
line), and 3.24 dB for state l|54[l (dotted line). The curves are 
plotted for r\ = 0.25. 



We plot the behavior of the fidelity and probability of 
generation of the target states l|51() - l|54|l as a function 
of the squeezing (Fig. 6), beam-splitter transmittancc 
(Fig. 7), and photodetector efficiency (Fig. 8). As in the 
preceding section, we observe that the fidelity of the gen- 
eration for any state gets arbitrarily close to one as T ap- 
proaches unity, as shown in Fig. 7. We also find that the 
fidelity is very robust against small detector efficiency rj, 
as can be seen in Fig. 8. On the other hand, the prepara- 
tion probability decreases with a growing T and decreas- 
ing 77, as predicted by the equation P oc (1 — T) 2 r] 2 . 

All these features are very similar to those found in the 
preceding section, where we considered only states gen- 
erated with one photon subtraction. Let us now stress 
some new features. First, we note here the existence of a 
clear optimal squeezing, giving the maximum fidelity for 
each of the four studied states, see Fig. 6(a). Observing 
that the optimal squeezing has a higher value [from 2.4 
dB for state l|53|l to 4 dB for state (|52f> ] than those en- 
countered in the case of one photon subtraction, we can 
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FIG. 8: (a) Fidelity between the generated state and the tar- 
get state and (b) probability of successful generation as a func- 
tion of r/ for the four target states (I51l > — (1541) . The curves are 
plotted considering the optimal squeezing si n for each state, 
namely 3.54 dB for state l|51|l (sòlid line), 4.02 dB for state 
(1521 (dashed line), 2.43 dB for state ||53|) (dot-dashed line), 
and 3.24 dB for state l|54[l (dotted line). The curves are plot- 
ted for T = 0.95. 



expect an increasing value of the optimal squeezing for 
an increasing number of Fock states in the target super- 
position. It can be checked that the value of this optimal 
squeezing tends to zero when T tcnds to 100%. 

Another interesting fact is the existence of very dif- 
ferent vàlues of the optimum fidelity for different target 
states for a fixed T = 0.95 and T) = 0.25, as shown in 
Fig. 6(a). For example, the squeezed two-photon state 
S(s)\2) is much more difficult to generate than the other 
three states ()52")) -(l54Tl. For the state S(s)\2), a trans- 
mittance of T > 0.99 is necessary to reach a fidelity of 
F > 0.95, rcsulting in a very low probability of gener- 
ation. This would make the experimental generation of 
5(s)|2) with a good fidelity very challenging. In con- 
trast, the squeezed balanced superposition state 1 (51)1 can 
be generated with a high fidelity F > 0.90 even with a 
transmittance T ps 0.90. 

Another surprising fact arises when /3i ^ p2 ■ Then, the 
equations H50[) give two distinct sets of a^s generating the 
same arbitrary target state, the second set being obtaincd 
by making the exchange (3\ «-> /?2 ■ Considering the pure- 
state description and T — > 1, the two alternative choices 
of displacements become strictly equivalent. In contrast, 
when considering the realistic model with T < 1, these 
two solutions for the same target state do not have the 
exact same behavior. As we can see in Fig. 9(a), one of 
the two solutions is indeed more robust to decreasing T. 
However, the two solutions are rather similar as far as the 
probability of state generation is concerned, as shown in 
Fig. 9(b). 



IV. CONCLUSIONS 

In this paper, we have shown that an arbitrary single- 
mode state of light can be engineered starting from a 
squeezed vacuum and applying a sequence of displace- 



F 




0.8 0.85 0.9 0.95 1 0.8 0.85 0.9 0.95 1 



FIG. 9: (a) Fidelity between the generated state and the 
target state and (b) probability of successful generation as a 
function of T for the two target states 152t and (1541) . The 
curves correspond to the two alternatives choices of the dis- 
placements a\ and Q2 when considering the optimal squeezing 
Sin for each state, namely 4.02 dB for state 15211 (dotted line, 
dot-dashed line), and 3.24 dB for state l|54[l (dashed line, sòlid 
line). The curves are plotted for r\ = 0.25. 

ments and single-photon subtractions. More precisely, 
the setup based on N photon subtractions can be used 
to generate a squeezed superposition of the N + 1 first 
Fock states. If this remaining squeezing is not desired, an 
extra squeezer can be added at the output of the setup 
in order to allow the generation of any (non-squeezed) 
superposition of the N + 1 first Fock states. 

We have shown that the desired target state can be 
successfully produced with arbitrarily high fidelity using 
a reasonably low squeezing 3dB) if the transmittance 
T of the beam splitter used for photon subtraction is 
sufficiently close to unity (e.g. T ~ 95%). This holds 
even when inefficient photodetectors with single-photon 
sensitivity but no single-photon resolution are employed, 
such as the commonly used avalanche photodiodes. We 
have studied the dependence of the achievable fidelity 
on the detection efnciency 77, and have found that the 
scheme is very robust in the sense that the fidelity does 
almost not depend on 77. Fidelities above 95% could be 
reached even with 77 of the order of a few percent if T 
is high enough. However, low r\ and high T drastically 
reduce the preparation probability, so that a compromise 
has to be made when determining T . 

Since our proposal does not require single-photon 
sources and can operate with low-efficiency photodetec- 
tors, its experimental implementation should be much 
easier than for the previous proposals, in particular the 
one based on repeated photon additions [Ï6J . The recent 
demonstration of single-photon subtraction from a single- 
mode squeezed vacuum 30] provides a strong evidence 
for the practical feasibility of our scheme. In this exper- 
iment, a rather low overall detection efficiency (rj w 1%) 
was reported, due to the spectral and spatial filtering of 
the beam measured by the photodetector. This is neces- 
sary because the OPA emits squeezed vacuum in several 
modes and the mode A is selected in the experiment by 
the strong local oscillator pulse in balanced homodyne 
detector. Despite the filtering, the single-photon detec- 
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tor PD can be sometimes triggered by photons coming 
from other modes, so this would probably be the main 
source of imperfections in our scheme. 

As a conclusion, we may reasonably assert that the 
generation of arbitrary squeezed superpositions of |0) and 
|1) should be experimentally achievable using our scheme 
with the present technology, while some improvement of 
the filtering and detection efnciency would probably be 
needed in order to extend the scheme to the preparation 
of squeezed superpositions of |0), |1), and |2). 
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